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I. Introduction 



The study of phase transitions in continuous spin lattice models has a long history. An 
important prototypical example of a random model in this class is the continuous spin random 
field model, where ferromagnetically coupled real valued spins fluctuate in randomly modulated 
local double- well potentials. 

In the present paper we study this model for weak disorder in dimensions d > 3 proving 
ferromagnetic ordering. Our aim is more generally to describe an expansion method mapping 
multiple-well continuous spin models to discrete spin models with exponentially decaying in- 
teractions by means of a single-site coarse-graining. Then we make use of information about 
the latter ones. This transformation can be regarded as an example of a useful (and moreover 
non-pathological) single-site 'renormalization group' transformation. While it is already inter- 
esting in a translation-invariant situation, it is particularly useful for non-translational invariant 
systems since it allows to 'factorize' the degrees of freedom provided by the fluctuations of the 
spins around their local minima. 

It is ten years now that the existence of ferromagnetic ordering for small disorder at small 
temperatures was proved for the ferromagnetic random field Ising-model (with spins a x taking 
values in {—1, 1}) by Bricmont-Kupiainen [BK1], answering a question that had been open for 
long in the theoretical physics community. The 'converse', namely the a.s. uniqueness of the 
Gibbs-measure in d = 2 was proved later by Aizenman and Wehr [AW]. For an overview about 
the random field model from the perspective of theoretical physics, see e.g. [Na]. Given the 
popularity of continuous spin models it is however certainly desirable to have a transparent 
method that is able to treat the additional degrees of freedom present in such a model. 

Bricmont and Kupiainen introduced in [BK1] the conceptually beautiful method of the 
renormalization group [RG] to the rigorous analysis of the low temperature behavior of a dis- 
ordered system, that turned out to be very powerful in this situation although there is no 
scale-invariance in the problem. The heuristic idea is: map the initial spin-system onto a coarse- 
grained one that appears to be at lower temperature and smaller disorder. Then iterate this 
transformation. This idea has to be implemented in a suitable representation of contours (that 
are the natural variables at low temperatures.) (For a pedagogical presentation of such a RG 
in application to the proof of stability of solid-on-solid interfaces in disordered media, see also 
[BoK], [K].) An alternative treatment of disordered lattice systems with finite local spin-space 
was sketched by Zahradm'k [Z2], however also using some iterated coarse graining. 

It is clear that also in the more difficult situation of continuous spins, spatial renormalization 
will be needed. However, continuous spins being more 'flexible' than Ising spins make it difficult 
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to cut the analysis in local pieces. It is then to be expected that the difficulties to control the 
locality of a suitably defined renormalization group transformation acting directly on continuous 
spins in a rigorous way would blow up tremendously compared with the discrete spin case of 
[BK1]. (The amount of technical work needed in their proof is already not small!) For an example 
of a rigorous construction of a RG-group for a continuous spin-lattice system, see [Bal], [Ba2] 
for the ordered Heisenberg-Ferromagnet. (This might give some idea of the complexities of such 
a method.) 

Indeed, despite the conceptual beauty, technical difficulties have kept the number of rigorous 
applications of the RG to low-temperature disordered lattice spin systems limited. Moreover, 
usually a lot of technical work has to be repeated when extending such a method to a more 
complex situation, while it would be desirable to make use of older results in a more transparent 
way. 

We will therefore describe a different and more effective way to the continuous spin problem: 
1) Construct a single-site 'RG'-transformation that maps the continuous model to a discrete one. 
Obtain bounds on the first in terms of the latter one. In our specific <fi 4 double-well situation this 
transformation is just a suitable stochastic mapping to the sign-field. 2) Apply the RG group to 
the discrete model. As we will show, the discrete (Ising-) model in our case has a representation 
as a contour model whose form is invariant under the discrete-spin RG that was constructed in 
[BK1]. So we need not repeat the RG analysis for this part but can apply their results, avoiding 
work that has already been done. 

In the last years there has been an ongoing discussion about the phenomenon of RG patholo- 
gies. It was first observed by Griffith, Pearce, Israel (and extended in various ways by van Enter, 
Fernandez, Sokal [EFS]) that even very 'innocent' transformations like taking marginals on a 
sub-lattice of the original lattice can map a Gibbs-measure of a lattice spin system to an image 
measure that need not be a Gibbs-measure for any absolutely summable Hamiltonian. (See 
[EFS] for a clear presentation and more information about what pathologies can and can not 
occur, see also the references given therein.) On the other hand, as a reaction to this, there has 
been the 'Gibbsian restoration program' initiated by the late Dobrushin [Do2] whose aim it is 
to exhibit sets of 'bad configurations' of measure zero (w.r.t. the renormalized measure) outside 
of which a 'renormalized' Hamiltonian with nicely decaying interactions can be defined. This 
program has been carried out in [BKL] for a special case (again using RG based on [BK1]). 

Since we will be dealing with contour representations of finite volume measures that provide 
uniform bounds on the initial spin system we do not have to worry about non-Gibbsianness vs. 
Gibbsianness to get our results. Nevertheless, to put our work in perspective with the mentioned 
discussion, we will in fact construct a uniformly convergent 'renormalized Hamiltonian' for the 
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measure on the sign-field, for all configurations. In other words, there are no pathologies in our 
single-site coarse graining and the situation is as nice and simple as it can be. 



Let us introduce our model and state our main results. We are interested in the analysis 
of the Gibbs measures on the state space = IR Z of the continuous spin model given by the 
Hamiltonians in finite volume A 

£™ A (m A ) 

= | ^2 ( m x - m y f + | ^ {m x -m v f + ^V(m x ) -^r] x m x (1-1) 

{",»)CA rtA; 9 e3A x £\ igA 

d(x,y) = l d(x,y) = l 

for a configuration m\ G Oa = IR A with boundary condition fhdh- Here we write dA = {x G 
A c ;3y G A : d(x,y) = 1} for the outer boundary of a set A where d(x,y) = \\x — y\\\ is the 
1-norm on IR d . q > will be small. Given its history and its popularity we will consider mainly 
the example of the well-known double-well (/> 4 -theory. As we will see during the course of the 
proof, there is however nothing special about this choice. We use the normalization where the 
minimizers are ±m*, the curvature in the minima is 1, and the value of the potential in the 
minima is zero and write 



V(m x , 



{m 2 - (m*) 2 ) 2 (L2) 



where the parameter m* > will be large. We consider i.i.d. random fields (i]x) xe z d ^ na ^ satisfy 

(i) r] x and —r) x have the same distribution 

(ii) IP[Vx >t]<e~^ 
(hi) \Vx\ <5 

where a 2 > is sufficiently small. The assumption (iii) of having uniform bounds is not essential 
for the problem of stability of the phases but made to avoid uninteresting problems with our 
transformation and keep things as transparent as possible. 

The finite volume Gibbs-measures [i™ 0A,VA are then defined as usual through the expecta- 
tions 



Ma 



(/) = ^rLl / dm A f(m A ,m Ac )e- E ^ VA ^ (1.3) 

£ A J IR A 

for any bounded continuous / on f2 with the partition function 

Z m 0A ,„A = f ,j fflAe -< 8A " A (».A) ( L4) 
J]R A 

We look in particular at the measures with boundary condition m x = +m* (for all x G Z d ) in 
the positive minimum of the potential, for which we write /i^ m ' VA ■ 
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To prove the existence of a phase transition we will show that, for a suitable range of 
parameters, with large probability w.r.t. the disorder, the Gibbs-expectation of finding the field 
left to the positive well is very small. Indeed, we have as the main result 

Theorem 1: Let d > 3 and assume the conditions (i) , (ii) , (Hi) with a 2 small enough. Then, 
for any (arbitrarily small) 7 > 0, there exist qo > (small enough), So, Si > (small enough), 
tq (large enough) such that, whenever 5 < So, q(m*) 2 > r and q(m*)i < Si we have that 



IP 



hmsup/i A 



m* ' 




const 


m x < ~Y 


>7 


<e - 2 



(1.5) 



for an increasing sequence of cubes An- 



Remark: Note that the quantity q(m*) 2 gives the order of magnitude of the minimal energetic 
contribution of a nearest neighbor pair of spins with opposite signs to the Hamiltonian (1.1); 
it will play the role of a (low temperature) Peierls constant. Smallness of q (to be compared 
with the curvature unity in the minima of the potential) is needed to ensure a fast decay of 
correlations of the thermal fluctuations around the minimizer in a given domain. The stronger 
conditions on the smallness, q < const (m*)~§, however is needed in our approach to ensure the 
positivity and smallness of certain anharmonic corrections. 

Let us now define the transition kernel T x ( ■ | ■ ) from IR to {—1, 1} we use and explain why 
we do it. Put, for a continuous spin m x 6 IR, and an Ising spin a x £ {—1, 1} 



T x (o x m x ^j :=-(! + a x tanh (am*m x )) 



(1.6) 



where a > 1, close to 1, will have to be chosen later to our convenience. In other words, the 
probability that a continuous spin m x gets mapped to its sign is given by ^ (1 + tanh (am* \m x \)) 
which converges to one for large to*. The above kernel defines a joint probability distribution 



H™ AC,VA (dmA)T (c?cta|toa) on IR A x {—1, 1} A whose non-normalized density is given by 



II T x{°~x\m x ) 



(1.7) 



Its marginal on the Ising-spins a\ 

(T ^t 0A ' VA )) (da A ) := ^ aA " A (dm A )T (da A \m A ) 
will be the main object of our study. 



(1.8) 



To prove the existence of a phase transition stated in Theorem 1 we will have to deal only 
with finite volume contour representations of (1.8), as given in Proposition 5.1. Nevertheless, 
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it is perhaps most instructive to present the following infinite volume result in the Hamiltonian 
formulation to explain the nature of the transformation. 



Theorem 2: Assume the hypothesis of Theorem 1 and let r/ be any fixed realization of the 



disorder. Suppose that fj, v is a continuous spin Gibbs-measure obtained as a weak limit of fij 10 



along a sequence of cubes A for some boundary condition fh G {— m* , m*} z . Then, for a suitable 
choice of the parameter a > 1 ( close to 1) in the kernel T the following is true. 

The measure T (fi 11 ) on {—1,1} Z is a Gibbs measure for the absolutely summable Ising- 
Hamiltonian 

Hl mn (a) 



I sing 

, A N-l * / A N-l ! / 

Vc 

x,y x C:\C\>2 



a 2 (m*) 2 v-^. » \ — x * / » \ — i ^ / n. (1-9) 
o l^( a - V A z*) x , y <7xO-y- am ^{a- qA Z d) xy r) x o x - ^ ^c(^c;^"' 



where A Z d is the lattice Laplacian in the infinite volume, i.e. A Z d. x ^ y = 1 iff x,y £ V are 
nearest neighbors, A Z d. Xty = —2d iff x = y and A Z d. XtV = else. 

The many-body potentials are symmetric under joint flips of spins and random- fields, &c(°~Ci r lc) 
^c(—o'Ci— r ]c), translation-invariant under joint lattice- shifts. They obey the uniform 
bound 

\$c(o-c,Vc)\<e-^ cl (1.10) 

with a positive constant 7. 



Remark 1: As in Theorem 1, 7 can be made arbitrarily small by choosing q , So, Si small and 
To large. More information about estimates on the value of 7 and 7 can in principle be deduced 
from the proofs. 

Remark 2: By imposing the smallness of S we exclude pathologies due to exceptional real- 
izations of the disorder variable 77 ('Griffiths singularities') in the transformation T. (We stress 
that this does not simplify the physical problem of the study of the low-temperature phases 
which is related to the study of the formation of large contours.) Starting from the joint dis- 
tribution (1.7) it is natural to consider the distribution of continuous spins conditional on the 
Ising spins; here the Ising spins a x will play the role of a second sort of external fields. Then, 
as it was explained in [BKL], possible pathologies in the transformation T would be analogous 
to Griffiths-singularities created by pathological Ising configurations. In this sense, Theorem II 
states that there are neither Griffiths singularities of the first type (w.r.t. n) nor the second type 
(w.r.t a). The treatment of unbounded random fields would necessitate the analysis of so-called 
'bad regions' in space (where the realizations of the random fields are anamolously large). This 
should be possible but would however obscure the nature of the transformation T. 
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Let us now motivate the form of T x and comment on the structure of the Hamiltonian. 
Introducing quadratic potentials, centered at ±m*, 

Q a *{m x ):=^{m x -a x m*) 2 + b (1.11) 

with b > (close to zero) to be chosen later, we can rewrite the transition kernel in the form 

T x (a x \m x ) = - q*.^.) (1-12) 

The crucial point is that the joint density (1.7) contains a product over x over the quantities 

e- yK) J , ,(^|m I )=e- Q " W (l + »K)) (1.13) 
where, using (1.12), we can write the remainder in the form 



-V(m x ) 



Now, if the initial potential V(m x ) is sufficiently Gaussian around its minima and the quadratic 
potential is suitably chosen, w(m x ) should be small in some sense. If w(m x ) were even 
zero, we would be left with OA-dependent Gaussian integrals that can be readily carried out. 
They lead to the first two terms in the Ising-Hamiltonian (1.9), containing only pair-interactions. 
This can be understood by a formal computation. The modification of the measure for 'small' 
w(m x ) then gives rise indeed to exponentially decaying many-body interactions, as one could 
naively hope for. 

Expanding n^eA (1 + w ( m x)) then leads in principle to an expansion around a Gaussian 
field. 1 One problem with this direct treatment is however that resulting contour activities will in 
general be nonnegative only if w(m x ) > for all m x . But note that the latter can only be true 
for the narrow class of potentials such that V(m x ) < Constm x for large \m x \. Thus, w(m x ) will 
have to become negative for some m x e.g. for V compact support or in the ^-theory. While it 
is not necessary to have positive contour activities for some applications (see [BChF],[Z3]) it is 
crucial for the random model: A RG, as devised in [BK1], needs non-negative contour weights. 2 



1 The author is grateful to M. Zahradmk for pointing out the idea to decompose e V{nn x ) into a sum of two 
Gaussians and a remainder term that should be expanded. However, contrary to [Z3] we write the remainder in a 

multiplicative form which allows for the transition kernel interpretation. 

2 

Vaguely speaking, the method keeps lower bounds on the energies of all configurations, but also upper bounds 
on the energies of some configurations (that are candidates for the true groundstates). This can be seen nicely in the 
groundstate-analysis of the models treated in [BoK]. To do an analogue of this for finite temperatures, non- negative 
(probabilistic) contour weights are necessary in this framework. 
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We are able to solve this problem and define positive effective anharmonic weights by a suitable 
resummation and careful choice of the parameters a, b of the quadratic potential Q ax ; these will 
be kept fixed. This choice is the only point of the proof that has to be adapted to the specific 
form of the initial potential V. Later the positivity of weights will also be used for the control 
of the original measure in terms of the Ising-measure (see Proposition 5.2). 

In Chapter II it is shown how non-negative effective anharmonic weights obeying suitable 
Peierls bounds can be defined. Chapter III finishes the control of the anharmonicity around the 
Ising model arising from the purely Gaussian theory (i.e. w(m x ) = 0) in terms of a uniformly 
convergent expansion. Chapter IV treats the simple but instructive case of the Ising field without 
the presence of anharmonicity, showing the emergence of (generalized) Peierls bounds on Ising 
contours. In Chapter V we obtain our final contour model for the full theory and prove Theorem 
1 and Theorem 2. The Appendix collects some facts about Gaussian random fields and random 
walk expansions we employ. 

Acknowledgments: 

The author thanks A.Bovier and M.Zahradm'k for interesting discussions and suggestions. 
This work was supported by the DFG, Schwerpunkt 'Stochastische Systeme hoher Komplexitat'. 

II. Anharmonic contours with positive weights 

We will explain in this Chapter how (preliminary) 'anharmonic contours' with 'anharmonic 
weights' that are non-negative and obey a Peierls estimate can be constructed. We start with 
a combinatorial Lemma 2.1. and a suitable organization of the order of Gaussian integrations 
appearing to derive algebraically the representation of Lemma 2.3. We will make no specific 
assumptions about the potential at this point that should however be thought to be symmetric 
'deep' double-well. Our later treatment is valid once we have the properties of 'positivity' and 
'uniform Peierls condition of anharmonic weights' that are introduced in (2.19) and (2.20). These 
are then verified for the (/> 4 -theory in an isolated part of the proof that can be adapted to specific 
cases of interest. 

We will have to deal with the interplay of three different fields: continuous spins m x (to 
be integrated out), Ising spins a x and (fixed) random fields r] x , subjected to various boundary 
conditions in various volumes. In some sense, the general theme of the expansions to come is: 
keep track of the locality of the interaction of these fields in the right way. For the sake of clarity 
we found it more appropriate in this context to keep a notation that indicates the dependence 
on these quantities in an explicit way in favor of a more space-saving one. 
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Now, since we are interested here in a contour-representation of the image measure T ^ / u™ SA ' 7 ' A ^ 
under the stochastic transformation (1.6), let us look at the non- normalized weights on Ising- 
spins given by 

Zf BA ' VA (a A ) := / dm A e- E Z dA ' VA ^ TT T x (a x \m x ) {2A) 



' ffiA x€A 

so that we get the desired Ising-probabilities dividing by Z^'^ = J2a A e{-i,i} A Zf aA ' VA (a A ). 

To describe our expansions conveniently let us define the following quadratic continuous- 
spin Hamiltonians, that are made to collect the quadratic terms that arise from the use of (1.13) 
to the above integral. We write, for finite volume VC/Z d , 

H rn av ,r, v ,a v 

= | ^2 ( mx ~ m y} 2 + I ( mx ~ ™y^ 2 + \ Z~2 ( mx ~ m * a x f ~ z^j Vx 771 * ( 2 - 2 ) 

{x,y}CV xGV ]y G8V xGV igV 

d(x,y) = l d(x,y) = l 

Here and throughout the paper we always write dG for the outer boundary inside A, i.e. dG = 
{x G A U B c ; d(x, G) = 1}. The notion 'nearest neighbor' is always meant in the usual sense of 
the 1-norm. The fixed Ising-spin ay € {—1, 1} V thus signifies the choice of the well at each site. 
From the point of view of the continuous fields it is just another parameter. 

With this definition we can write the non-normalized Ising- weights (2.1) in the form 



Z™ M ^((7 A ) = e" b l A l f dm A e-< 8A,,A " A ( mA ' TT (l + w(m x )) 

J ^ A reA 



(2.3) 



If the w(m x ) were identically zero, we would be left with purely Gaussian integrals over Ising-spin 
dependent quadratic expressions. This Gaussian integration can be carried out and yields 

/ dm A e-<° A ^ A ^ = C A x e- in( ^ A K 8 ^^) {2A) 
Jm A 

with a constant C A that does not depend on a A (and r] A ). The latter fact is clear since a A 
(and 7] A ) only couple as linear terms ('magnetic fields') to m A while they do not influence the 
quadratic terms. Note the pleasant fact that no spacial decomposition of the Gaussian integral 
is needed here and no complicated boundary terms arise. 

Now the minimum of the continuous-spin Hamiltonian in the expression on the r.h.s. of (2.4) 
provides weights for an effective random field Ising model for the spins a a', its (infinite volume) 
Hamiltonian is given by the first two terms in (1.9). The treatment of this model is much simpler 
than that of the full model; all this will be postponed to Chapter IV. There it is discussed in 
detail how this model can be transformed into a disordered contour model by a mixed low- 
and high-temperature expansion. However, since this model provides the main part of the final 
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contour model that is responsible for the ferromagnetic phase transition some readers might 
want to take a look to Chapter IV to understand the form of our final contour-representation in 
a simpler situation. 

Our present aim now is however to show how the anharmonic perturbation induced by the 
u>-terms can be treated as a positive-weight perturbation of the purely Gaussian model. 

Let U = U + U {—U + )dIR, where U + is a suitable 'small' neighborhood of the positive 
minimizer of the potential m* that will be determined later and that will depend on the specific 
form of the potential. The first key step to define non-negative activities is now to use the 
following combinatorial identity on the set U = {x € A; m x £ U}. 

Lemma 2.1: Let k<Z/Z d be finite and connected. For any setUdk we can write the polynomial 
nLeA (1 + w x) i- n the |A| variables (w x ) xeA in the form 

=i + e n n w 

G:0/GcA G* x€dd 

conn, cp of a 



n (WW + w x ) - n 



^x<fU 



(2.5) 



The proof is given at the end of this chapter. Application of Lemma 2.1. gives us the 
expansion 



Jm A 



'1R A 



G:^GGA JlRA (2.6) 



Gi x€dGi 

conn.cp of a 



H 0-m x eu + w{m x )) - \ \ 

.x€Gi x€Gi 



Note that the expression under the integral factorizes over connected components of G := GUdG. 

To introduce the anharmonic (preliminary) weights we need a little preparation. To avoid 
unnecessary complications in the expansions it is important to organize the Gaussian integral in 
the following conceptually simple but useful way: We decompose the nonnormalized Gaussian 
expectation over the terms in the last line into an outer integral over rridG and a 'conditional 
integral' over m^\Q G given mgo. The latter integral factorizes of course over connected compo- 
nents of A\dG; in particular the integrals over A\G and G become conditionally independent. 
W.r.t. this decomposition they appear in a symmetric way. 
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To write down the explicit formulae we need for that we introduce 

Some notation: The V x ^-matrix Ay is the lattice Laplacian with Dirichlet boundary- 
conditions on VcA, i.e. Av. x ,y = 1 iff x,y £ V are nearest neighbors, Av- X . y = —2d iff 
x = y £ V and Av ]X , y = else. II y is the projection operator onto Oy (in short: onto V), i.e. 
^-v-,x,y = l x = y ev- We also use the redundant but intuitive notations m A \v = Iiv m \ = m v f° r 
the same thing. ly is the vector in IR A given by ly, x = Ixev- For disjoint Vi, V2CA we write 
dv 1 ,v 2 for the matrix with entries dv lt v 2 - x , y = 1 iff x £ V\, y £ V 2 are nearest neighbors and 
dv 1 ,v 2 ;x, y = else. We write Ry ■= (c — Ay) -1 for the corresponding resolvent in the volume 
V. Here and later we put c = -. 

For the sake of clarity we keep (at least for now) the dependence of all quantities on 
continuous spin-boundary conditions, random fields, Ising-spins, as superscripts. Then we have 

Lemma 2.2: For any subset GcA the random quadratic Hamiltonians (2.2) have the decom- 
position 

A Tjrri3A,VA,0'A I \ 1 A IT™M > m 8G >VA\dG i CT A\ 8G / \ . • r TTlflg A ,77A ,<J A I I \ (~ 

= AH dSX (rn dG ) + AH A ^ dG {m A \ dG ) + mf H A SA ' (m A ) 



Here the 'fluctuation- Hamiltonians' are given by 

AH rnaA,r,A,*A 

= \< (m dG - mt B *™'° A \ dG ) , (n 9G (a - qA^ 1 II 9G ) {m dG - mf^^\ dG ) >*? 

(2.8) 

and the 'conditional fluctuation- Hamiltonian' (i.e. conditional on rriQ G ) 

A\dG {rn A \ dG ) 
= 2 \ mA \ dG ~ m A\dG ) ' ( a ~ Q A A\dG) [m A \ dG - m A ^ dG J > A \ G 

As centerings are occuring: the 'global minimizer' 



m 



mgA,riA,<7A 
A 



(2.9) 

R A ^cm*a A + -^- + d^dArhdr^j ( 2 -10) 



and the 'conditional minimizer' 

m dA ,m G,riA\aG^A\aG 
m A\dG 



^A\9G [cm (Ja\9g H 1- o A \ dGtdG m dG + o A \ dGtdA m dA 



(2.11) 



The proof is a consequence of Appendix Lemma A.l(iii) which is just a statement about 
symmetric positive definite matrices. Lemma 2.2 can be seen as an explicit expression of the 
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compatibility property for the Gaussian local specifications denned thru the Hamiltonian (2.7) 
in the volumes A\<9GcA. Indeed, the Gaussian measure defined with the quadratic form (2.8) 
describes the distribution on A projected onto dG. (Since we will use this formula later for 
subsets of A it is convenient to make the A explicit at this point, too.) The Gaussian measure 
on A\dG defined with (2.9) is the conditional measure given maG- 

We like to stress the following decoupling properties of the conditional expressions. Equation 
(2.11) for the conditional minimizer decouples over connected components of A\dG since the 
resolvent Ra\qg is J us t the direct sum of the Rv t 's. So we have that 



m A\dG \ Vi = R Vi \ cm a Vi + — + Ov^dVi-maVi + ov^dKmah I 



--:m Vi 



(2.12) 



is a function depending only on what is appearing as superscripts, namely random fields and 
Ising-spins inside Vi and continuous-spin boundary condition on dVi. (The dependence on the 
global boundary condition toqa is of course only thru rh x for d(x, Gi) = 1. We don't make this 
explicit in the notation.) 

Also, the conditional fluctuation-Hamiltonian on A\dG decomposes into a sum over con- 
nected components of its support A\dG: 



a\og { m A\dG) = 2^ AH Vi ( m vJ where 

i 

AH™; A ' maVi ' Vv *' aVi (m Vi ) (2-13) 

1 . / rh gA ,m gv ,rj Vi ,o- Vi \ i a \ ( mg A ,mg V ,-qv ,(T Vi \ 

= g < [ m v t ~ m v t J > (a - Q A vJ [m Vi - m v . ) >v * 



Putting together the connected components of A\G we can thus write 

AH ^ G rng GWW g a 

= AH^ d A G mda ' V ^' a ^ ( TOax5 ) + £ — ^ (ma,) (2-14) 

conn.cp of a 



So, the sum over G"s in (2.6) can be written as 
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dm^e 



G:0/GCA 



x n n im ^ u 

conn.cp of g 

G:0/GcA 



— A_H' maA '* 7A '' JA 

dmQce SG - A 



'" A ' CTA (-a G ) TT i jr 
11 l m x £U 

x£dG 



d m A\G e AXG 



( m A\c) 



n 



_m gA ,mg G ,77 G ,cr G 



(m Gj ) 



conn.cp of g 



11 O-mxtu + w{m x )) - \\ l mx ?u 

.x€Gi x€Gi 



(2.15) 



Now we note the pleasant fact that the Gaussian integral over A\G is independent of all of the 
superindexed quantities (since they appear only in the shift of the quadratic form), so that it 
can be pulled out of the mac-integral. It gives 

-™e A ' m 8G' ,J A\G' C 'A\G/ \ 
\ m A\G/ 



/ 



-AH .... 

dm A\G e MG 



( 2.)^(det(a-,A A ^))" 5 (2.16) 



Let us look at the last line now. Conditional on mgc we define anharmonic activities by the 
formula 



T m gA ,m gG 



> r IO i ,<TG i ._ J 



dm G .e 



-AH. 



{ m Oi) 



x€Gi x€Gi 

(2.17) 

We write j'^ 6A ' maG ' riG '' TG = i f or (7 = 0. So we have obtained the following representation for 
the non-normalized Ising-weights 

Lemma 2.3: With the above notations we have 

e6 |A|^a A(aA)=e -^<— K) 

E (27r)^(det(a-gA A ^))" i 



G:0CGCA 
, -AH^^^imsa) TT -, 

dm aG e sg,a v bg, || l m;cer7 

xeaG 



(2.18) 



n 

conn, cp of a 



r™9A ,meGi j^G; ,<7Gi 

1 G l 



Let us pause for a minute and comment on what we have obtained. For the purely Gaussian 
model (i.e. the -w-terms are identically zero) the contributions for G ^ vanish. So the above 
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formula is a good starting point for the derivation of the signed-contour representation whose 
main contributions are provided by the minimum of the Gaussian Hamiltonians in the first line. 
The main other non-trivial ingredient are the preliminary anharmonic activities i^9A,m aa ,r] G ,a G ^ 
First of all, the whole construction makes only sense, if we are able to prove a suitable Peierls 
estimate for them, to be discussed soon. They factorize over connected components of the 
set G. The conditioning on mdG has allowed us to have them local in the sense that they depend 
only on random fields and Ising-spins inside Gi . Note that such a factorization does of course not 
hold for the remaining integral over dG (that would mean: over connected components of dG), 
as it is clear from (2.8). Indeed, the fields mgc fluctuate according to the covariance matrix in 
the total volume A. So to speak, their (stochastic) dependence is mediated by the Gaussian local 
specification defined with (2.8). Furthermore, the dependence of their mean-value in this local 
specification is (weakly) on all Ising-spins and random fields in A. Both kinds of dependence 
will have to be expanded later in Chapter III when the integral over dG is carried out. This will 
be done by enlarging the 'polymers' G and performing a high-temperature expansion. Finally, 
the determinants provide only trivial modifications of the weights that we will obtain; they can 
easily be handled by a random walk expansion. 

Let us stress the following nice feature of the above representation: 'Low-temperature 
contours' (see Chapter IV) will be created only by the global energy-minimum in the first line. 
Consequently there will be no complicated boundary terms for these 'low-temperature' terms 
(that could be easily produced by a careless expansion). 

Our further treatment of the expansion will be done under the assumption of the following 
two properties: 

Positivity of anharmonic weights: 

jrh 8A ,m gG ,ria,^a y q ^2 19) 

for all connected G, and m dA £ U dA , m dG £ U dG , n G £ [-5,5]°, a G £ {-1, 1} G . 
Uniform Peierls Condition for anharmonic weights: 

jm gA ,m gG ,riG,VG < £ |G| ^ 20) 

for all connected G, and rhdk £ U dA , mgc £ U dG , tjg £ [— S, 5] G , o~g € {—1, 1} G with e > 0. 



Rather than trying to be exhaustive in the description of potentials that satisfy these 
conditions we will use the rest of this Chapter to fix some properties that imply them and 
discuss in detail the explicit example of the 4 -theory in Lemma 2.6. This should however 
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indicate how the above two conditions can be achieved in concrete cases by suitable choices 
of the neighborhood U and the constants a and b occuring in the quadratic potential. The 
expansion will be continued in Chapter III. 

Let us start by fixing the following almost trivial one-site criterion. It makes sense if we are 
assuming the nearest neighbor coupling q to be small. 

Lemma 2.4: Suppose that w{m x ) > for m x G U. 

(i) Assume that we have uniformly for all choices of superindices 

dm x e 2 v - - ) w(m x )l mx€U 

(2.21) 

dm x e *y m * m * ) l m ^u 

Then we have the positivity (2.19). 

(ii) Assume that 

dm x e m * ) {w{m x )l m ^u + {I + w{m x ))l m ^u) < e (2.22) 

Then we have the uniform Peierls estimate (2.20) with the same e. 
Proof: Since we always have — 1 < w(m x ) < oo the assumption l mx€ uw(m x ) > implies that 

duiQe g (™g) || (l maiu + w (m x )) 

x€G (2 23) 

> / dmae-^G^^ JJ w{ m x )l mxeU > 

We reduce the estimation of the integrals to product integration by the pointwise estimate on 
the quadratic form 

a\\v G \\ 2 2 << v G , (a - qA§) v G > G < (a + 4dq)\\v G \\t (2.24) 

This gives 

/ dm G e \ m G ) ^ w (m x )l maeu 

■> x€G 

/ dm x e 2 y m * m * ) w{m x )l mx€U 



(2.25) 



and, on the other hand, 



7 _A H rh 8A<™3G- r lG-<*G ( ) T-r 

dm G e g (™g) || lm ^ v 



J 

/ dm x e 2 ^ m * ) l m ^u 

... r- J 



(2.26) 



x€G 
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This proves (i). 

The Peierls estimate (ii) follows from dropping the second product in the definition of / 
and using (2.24) to write 

Tm gA ,m gG ,r] G ,a G 
1 G 







< f d m G e- AH c aA ' maa ' va " a (™e) Yl ((l + w (m x ))l m ^u + w(m x )l mx€U ) 

^ x€G 
TT f a (m rn fh dA< m dG' r IG'< J G\' 2 

< I I / dm x e-*V"--"" ) ({l + w(m x )) + w(m x ) 



(2.27) 



Next we compute how big the nearest neighbor coupling q and size of the random fields 5 can 
be in order that any boundary condition in U yields a minimizer of the Gaussian Hamiltonian 
on G that is 'well inside' U. We have 

Lemma 2.5: Let < A± < A 2 and U + = [m* - A 2 ,m* + A 2 ], U = U + U (-U + ). Assume 
that q<^ ( K 2m * A + 1 A2 - l) and5 < zAl. Then we have that 

\ m fBA,m gG , VG ,a G _ m * ax \ < Ai (2.28) 
for all G, m aA e U dA , m aG G U dG , Vg G [-5,6]°, o G G {-1, 1} G . 

Proof: Note the linear dependence m fdA,m gG , VG ,a G = m m gA , mgG , VG =o,a G + (jl G ^ . Let us 
thus choose the condition for q s.t. 



rn 



m gA ,m gG , Tj G =0,<T G 



m*a x \ < ^ (2.29) 



This condition is in fact achieved for a one-point G = {x} and the boundary conditions having 
the 'wrong sign' with modulus m* + A 2 as we will formally see as follows. Let us assume that 
a x = — 1 and write this time for simplicity dG for the boundary in Z£ d (including possible sites 
in the outer boundary of A in /Z d ). Then we have, due to the positivity of the matrix elements 
of R G that 



m rh gA ,m gG , r] G =0,a G 



< -R G;x , x cm* + R G ; X , y cm* + (R G d G , dG l aG (m* + A 2 )) x ( 2 - 30 ) 

V€G\{x} 

We employ the equation R G (cl G + d G ^ G \ aG ) = l G to write the last line of (2.30) as 

m* - 2R G;XjX cm* + A 2 - A 2 (R G cl G ) x (2.31) 
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We note that Rg- x ,x is an increasing function in the sets G B x (which can be seen by the 
random walk representation, see Appendix (A. 8)). Further (Rg^g) x ls an increasing function in 
G. So the maximum over G of (2.31) is achieved for G = {x}. With R{ x y- tXtX = the value 
of (2.31) becomes -m* + (2m* + A 2 ) ^ which gives the upper bound 



m. 



dA,moa,VG=0,Cr G \{x} ,<Ta; = — 1 



+ m* < (2m* + A 2 )-^ (2.32) 

In the same way we obtain 

m x +m > -A 2 — — -T (2.33) 

c + 2d 

Equating of the r.h.s. with Ax/2 gives the r.h.s. of the condition on q stated in the hypothesis. 

For the estimate of the random field term note that < Rc;x,y < Rz d ;x, y an d Yl y ez d R z d ;x, y = 
- which give us 



yeG q 



5 \ -> „ $ A-i , 

<- E **'i*,v = a - T (2 ' 34) 



<> 



At this stage the treatment has to be made specific to the concrete potential and we spe- 
cialize to our example, the 4 -theory with potentials given by (1.2). The following Lemma 
summarizes how we can produce positivity and an arbritrarily small anharmonic Peierls con- 
stant. More specific information can be found in the proof. 

Lemma 2.6: For fixed e > we put 

U+ = [m* - (e m*)^ ,m* + (e m*)*] (2.35) 

Then we have 

(i) For any value of eo,m* , q, 5 there exists a choice of parameters a and b such that the an- 
harmonic weights obey the positivity (2.19). 

Furthermore there exist strictly positive constants a(m*, eo), b(m*, eo), qo(m*, eo), and 5o(m*, eo) 
such that the following is true. 

(ii) For all q < qo(m* , e ) and 5 < So(m* , eo) we have the Peierls estimate (2.20) with a constant 
e(eo,m*) that is independent ofq,5. 

(Hi) If 'eo is small enough this constant obeys the estimate e(e ,m*) < |g whenever m* > m^eo) 
is large enough. 
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The above constants can be chosen like 

, (2 + epW-g) 2 
a(m ,e ) = ~ 1 

(2.36) 

/ * x a(m*,e ) / i 2 X" 1 a(m*, e ) (e m*) 3 
g(m,e ) = (^20e 3 m 3 + gj ; 5 Q (m,e ) = 



20 



and 6(m*,e ) ~ e - co ™ s * m * 3 with m* | 00. 



Proof: We will take time to motivate our choices of the parameters that are made to ensure 
the validity of the assumptions of Lemma 2.4. Let us write the neighborhood U + in the form 
U + = [(1 — ei)m*, (1 + e\)m*\ and show why the choice of e\ given in (2.35) comes up. The 
zeroth requirement on a and b we have to meet is w(m x )l mxe u > 0. So, let us choose the 
Gaussian curvature a > 1 to be the smallest number s.t. we have, for all m x G U + , that the 
Gaussian centered around m* is dominated by the true potential i.e. 

e - « m *- m ' )2 < e -v(m B ) (2.37) 

with equality for m x = (1 + ei)m*. This amounts to a = ^ 2+ l^ , as in (2.36). Then we have on 
U + for the Gaussian centered around — m* 

e - <^+™*) 2 < e _ (^x) 2 ( 2 -n) 2 + i-(i+n) 2 m »2 e _ y(mx) (2 38) 

which gives us the estimate 

1 + w(m a; ) > e b 

on Any choice of e b bigger than the denominator thus ensures w(m x )l mxe u > 0. 



(2 + ei ) 2 (2- ei ) 2 + l-(l + ei ) 2 _»2 

1 + e s m 



(2.39) 



To have property (i) in Lemma 2.4. we have to choose e b even bigger. Obviously it is 
implied by 

inf WmJ > : N2 (2.40) 

But note that we always have 

\\mf aA ' maG ' VG ' aG \-m*\ <m max (m*,5,q,a) (2.41) 

with a constant rn max (m*, <5, q, a) that is finite for any fixed m*,S, q, a and that is estimated by 
Lemma 2.5. So the trivial choice 

e Hrn*,6, q ,a) . = ^ + - < 2 +n> 2 ( 2 -n) 2 + i-d + n) 2 m „ 

x fi + s,p '^r^V*" ) (2 ' 42> 

Y m-.\rh\<rh™^( m * ,S,q,a) Jdm x e 2 m ) \ mx& j) 
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gives some finite number and ensures the positivity of the anharmonic activities. This proves 

Let us now turn to quantitative estimates on the Peierls constant. To start with, the above 
definition of b is of course only useful if b will be small. Now, the r.h.s. of (2.42) is small 
whenever the centering of the Gaussian integrals is 'safe' inside U and the neighborhood U is 
big enough to carry most of the Gaussian integral. We apply Lemma 2.5. with A 2 = eim* and 
A 1 = jft. The hypotheses of the Lemma then give us the conditions q < qo and 8 < Sq with 



a /20 



So 



Then we have 
f dm x e 



2" I " L x — "t x 



) 1 



m x ^U 



j dm x e 



''X n^x 



rn gA ,m gG , VG ,v G 



) 



< 



a + 4dq 



ae\m 
20 



IP 



(2.43) 



|G| > 



10 



a 1-IP[\G\ > Va + 2dq 



9eim* 
10 



(2.44) 



This shows that b ~ e~ const ' (tim "^ tends to zero rapidly if eim* is getting large. 

Let us now see what Peierls constant we get according to Lemma 2.4 (ii). This will explain 
why the neighborhood U + should in fact be of the form (2.35). 



Our choice of U and a yields that we have, for all m x £ U + , that 

e -V(mx) + — 2 < e e i( m x-m ) 

This gives 1 + w(m x ) < e b +£i( m x- m ) . From this we have 

(lin,' I (1 + w(m x )) \ mx&J ax 

\ 2 „ 

) +ei(m x -m ) 



(2.45) 



< e b I dm r e 



f (m x - 



■mg A ,m ga , VG ,cr G 



(2.46) 



= e 



p «-2q \ m * 



m OG' r IG' c 'G . 



a — 2ei 



< e l 



2ir 



a — 2ei 



.g 100(a-2 ei ) 



and hence 



/ 



dm x e 



m gA ,mg G ,r] G ,cr G 



)'<! 



(1 + iy(m x )) l ma;e c/ 



< 2e b 



2vr 



_gl00(a-2 ei ) 



2tt 




iP 



|G| < 



9eim* 
10 



(2.47) 



a — 2ei 

Indeed, the l.h.s. is 0{e\m* 2 ) + O(ei) and thus imposes the condition that e\m* 2 be small! 
This estimate can essentially not be improved upon. It determines the dependence of the Peierls 
constant e on ei and m*. 
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Finally, the integrals over U c are much smaller: Indeed, for the bounded part of U c we 
estimate 



L 



< 



(l-ei)m* 



ma A ,m 



dm T e 



(l-ei)m* 



dm T e 



•vg>°g\ 

) (l + w{m x )) 

/•(l-ei)m „,'™9A> m 9G' , 'G>< J G N \ /•_ „,*> „, \ 



(2.48) 



We have for the last integral 



<-L 



(l-ei)m* 



dm x e 



-a (r 



ma. ,mc 



(1— ei)m* 



oei m* , * x (m x - m 



(2.49) 



The maximizer of the last exponent is m x = m* + 2ae 1 1 m which is outside the range of integration 
(due to our choice of the 10 before (2.43)) Estimating for simplicity the integral by the value of 
the integrand at (1 — e\)m* just gives 



/ (1 " ei)m ^rn, e -t(--- SA ' maG ^^)V + -(^))<m* e -(i-^)(-^) 2 
Jo 



(2.50) 



For the unbounded part of of U c where m > m*(l + ei) we have with our choice of a that 
1 + w(m x ) < 1. This gives us 



POO 



dm x e 2 



ma a ,mc 



'(l+ei)m* 




< W— IP 

a 



9eim* 
10 



(2.51) 



<^ g — const (eim*) 2 



Collecting the terms gives our final estimate on the Peierls constant 



e < 2e b 



2vr 



2vr 



glOO(a-2 £1 ) _ , / 

a — 2ei V a 




+ m * e -(!-^™ ) 2 +3 A /— IP 

a 



9eim* 
10 



(2.52) 



From here the lemma follows. 



We are still due the 
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Proof of Lemma 2.1: We expand ELeA C 1 + «>*) = 1 + E Ao :0/A n cA EUao Let 
-A(A )c(A\W)\A denote the maximal set amongst the sets Ac(A\U)\A that are connected to 
A . (We say that a set A is connected to a set A iff, for each point u in A, there exists a nearest 
neighbor path inside A U A that joins u and some point in A .) Equivalently, this A(A ) is the 
unique set AcA\A s.t. x U for all x € A and x G W for all x G <9(A U A). 

We collect terms according to the sets G = A U ^4(A ). Denoting by Gi the connected 
components of G and by Lj = A n Gj we have then 

jj (i + = i + yi n x *z u n n w - 

i6A A o :0#A o CAzeA(Ao) rEeO(A U J 4(Ao)) xeA 

= i+ e n En w« n wn 

G:0/GcA G i Lr.^LiCGi xeGi\Li x€dd 

conn.cp of g 

Adding and subtracting the term for = we have 

n n ^ = n n (2.54) 

Li-.tl^LiCGi x€Gi\Li xeLi xed x€Gi 

which proves the lemma. <C> 



(2.53) 



III. Control of Anharmonicity 

We start from the representation of Lemma 2.3 for the non-normalized Ising weights. We 
assume positivity and Peierls condition for the anharmonic (/-) weights as discussed in Chapter 
II and verified for the <f> -potential. Carrying out the last remaining continuous spin-integral we 
express the last line in (2.18) in terms of activities that are positive, obey a Peierls estimate and 
depend in a local way on the Ising-spin configuration o~\ and the realization of the random fields 
77a- We stress that all estimates that follow will be uniform in the Ising-spin configuration and 
the configuration of the random field. 

The result of this is 

Proposition 3.1: Assume that the anharmonic I -weights (2.17) satisfy the Positivity (2.19) 
and the uniform Peierls Condition (2.20) with a constant e. Suppose that e is sufficiently small, 
q is sufficiently small, a is of the order one, q(m*) 2 sufficiently large. Suppose that 5 < Const m* 
and \U\ < Const m* with constants of the order unity. 

Then, for any continuous- spin boundary condition itiqa £ U dA o,nd any realization of the 
random fields n\ G [— 5, 5] A , the non-normalized Ising weights (2.1) have the representation 

,iai l A l 1 - inf , ff" 1 SA>'7A.''A ( m l \ 

^A,*7A ((7A) = e -6|A| (27r) y (det(a-gA A ))-* e ™a H * ™ 

G:0CGCA 
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where the activity p appearing under the G-sum is non-negative and depends only on the indicated 
arguments, p factorizes over the connected components Gi of its support G, i.e. 

jTW (G; o- G ,7] G ) = YlfTWi (G i; a Gi , VGi ) (32) 

i 

and we have /f 8 sa g (G = 0; o- G ,7] G ) = 1. 
p has the 'infinite volume symmetries' of: 

(a) Invariance under joint flips of spins and random fields p (G; a G ,7] G ) = p(G;—a G ,—r] G ) if 
G does not touch the boundary (i.e. dg\G = 

(b) Invariance under lattice shifts p (G;o G ,rj G ) = p (G + i; a G+t ,VG+t) if G, G + tcA don't 
touch the boundary 

We have the uniform Peierls estimate 

r^{G;a G , VG )<e- a ^ (3.3) 

with a = const x min |log ^, log ^ ^ \ogm* ) | ■ 

Remark 1: Note that the first line of (3.1) gives the value for vanishing anharmonicity (i.e. 
w(m x ) = 0). 

Remark 2: For any fixed Ising-spin a a and realization of random fields rj\ the sum in the last 
line is the partition function of a non-translation invariant polymer model for polymers G. Note 
that there is no suppression of the activities p in the above bounds in terms of the Ising-spins. 
From the point of view of the polymers G the Ising spins and random fields play the similar role 
of describing an 'external disorder.' 

Proof of Proposition 3.1: To yield this representation we must treat the last line of (2.18). 
We can not carry out the mQ G -mkegYsl directly but need some further preparation that allows 
us to treat the 'long range' parts of the exponent by a high-temperature expansion. Depending 
on the parameters of the model (to be discussed below) we will then have to enlarge and glue 
together connected components of the support G . For any set GcA we write 

G r = {x G A; d(x, G) < r} (3.4) 

for the r-hull of G in A. Then we have, under the assumptions on the parameters as in Propo- 
sition 3.1. 

Lemma 3.2: There is a choice of ' r ~ Const such that the following is true. For each 

fixed subset GcA, continuous- spin boundary condition nig^ G U dA , fixed Ising- configuration 
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cja G {—1, 1} A and random fields t]a € [— 5, 5] A we can write 
' dmeae-^X^^o) ^ 



(2tt)^ Jde^n^a-gAGO-'lW) £ (g, G; a & , v G ) 



(3.5) 



G:GCA 
G'CG 



where the activity appearing under the G-sum depends only on the indicated arguments and obeys 
the uniform bounds 

< p™^o (g, G; ag, r^) < e~^\ (3.6) 

i/rti/i a = Const x min |log |,log \ (1^7^^ |- It factorizes over the connected components Gi 
of the set G, i.e. 



For G not touching the boundary (i.e. dg\G = $) p is invariant under joint flips of spins and 
random fields and lattice shifts. 

Remark: Later it will be convenient to have the determinant appearing on the r.h.s.; in 
fact it could also be absorbed in the activities under the G-sum. 

Proof of Lemma 3.2: Let us recall definition (2.8) of the 'fluctuation- Hamiltonian' 
(involving the global minimizer (2.10)) which gives the Hamiltonian of the projection onto dG 
of an Ising-spin and random-field dependent Gaussian field in A. Our first step is to decompose 
this projection from A onto dG into a 'low temperature-part' and a 'high temperature-part'. 
For fixed G we will consider definition (2.8) where A will be replaced by G r ; for r large enough 
the resulting term 'low-temperature'- term is close enough to the full expression, so that the rest 
can be treated by a high-temperature expansion. 

We write &b ■= {x G B;d(x,A) = 1} for the outer boundary in a set BdZ, d . Recall that, 
with this notation dA = d A A, so that d zd {G r ) = d dA {G r ) U d{G r ). 

Then the precise form of the decomposition we will use reads 

Lemma 3.3: With a suitable choice of r ~ Const lo g we have 

Mi) 

A#£?r A ' CTA (mac) 

= AH^Sr'° dAGr) ' VGr ' aGr (m dG )+ £ m T G M m 0G,a Gr , VGr ;C,a c , Vc ) ( 3 - 8 ) 

CCA 

CnSG#0;Cn(G T ') c #0 
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where the functions appearing under the C-sum depend only on the indicated arguments and 
obey the uniform bound 



\H^M m BG^Gr,VGr;C,a c ,Vc)\ < e" 5|C| (3-9) 

uniformly in rriQG € U dG and all other quantities for the C's occuring in the sum in (3.8). Here 
a = const log - . 

Remark: Note that the first part ('low temperature-part') decomposes of course over the 
connected components (G r )i of G r , i.e. 

AH dG,G- \ m dG) ~ 2^ /XH dGn(G-) z ,(G-) z {m gGn ( G r).) (3.10) 



Proof of Lemma 3.3: The l.h.s. and the first term on the r.h.s. of (3.8) differ in two 
places: The matrix and the centerings. We expand both differences using the random walk 
representation. 

The decomposition of the matrix into the matrix where A is replaced by G r and a remainder 
term can be written as 

(UaGRAUaG)- 1 

= (UecRGr-UeGy 1 - d dGA \ dG Tl(- -► ■■,C)d A \ dG , aG ( 3 - n ) 

CCA\8G 
Cn(G r ) c ^0,CnG 2 7i0 

where the A x A-matrix 1Z (• — ► ■ ; C) has non-zero entries only for x, y £ C that are given by 



oaths "v from x to u v ' 



paths -y from a; to y 
Range(7)- C 



(3.12) 



For the proof of this formula see the Appendix (A. 8) and (A. 13) ff. where also more details 
about the random walk expansion can be found. 

Simply from the decomposition of the resolvent R A = R G r + (Ra — Rc r ) an d the random 
walk representation for the second term follows the formula for the centerings 



,ma A ,r) A ,a A _ ^mg A ,T] a r ,a G r 

' (3-13) 



CCA 

cn(G'') c #0 



with 

m 



™eA,VGr,*Or . = ^ ^ ^ + ^ + d G r tdA m dA ^ (3.14) 
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and 'high-temperature' terms given by the matrix product 

m{C\ ac,Vc) = K (• -»■ ■ ; C) ^cm*a A + ^ + 5 c ,aA?Ti 9A ^ (3.15) 
From the bound on the resolvent (A. 12) we have uniformly 

\rh x (C;a c ,Vc)\ < Const (m* + 5) (l + — j (3.16) 

This quantity is in turn bounded by, say, ( 1 + 2§~ ) if we have that \C\ > r with r := 

Const - — l ? s m - t- . So we have r ~ Const |° s /"? \ for small o. 

Mi) q 

To write both type of summations over connected sets C in the same form we note that 



^ ddG,A\dG^(- — > - ; Cl) dA\dG,dG 

C 1 CA\3G 

c 1 n(C)° 7 ;0,c 1 nG 2 7 i0 ^ ^ 

^ ddG,A\dG^(- — ► • ; C 2 \dG) d A \ dG:dG lc 2 \dG 



C 2 CA 

c 2 neG^0;C2n(G r ) c 7;0 



which gives us the same range of summation for both sort of terms. The expansion then produces 
triple sums over connected sets C. Collecting terms according to the union of the occuring C's 
we obtain the desired decomposition with 

HdG,G- ( m dG , o-Q- , Vg- ;C,a c ,7] C ) 
q 

= -- < (m dG - m a dG r ) ,d dG)A \ dG K(- -> -;C\dG) d A \ dGydG l c \dG conn, ('mac ~ rrf dG ) > 
+ q < (m dG - mg G r ) ^UdcRcrUdcy 1 fh(C;o- G ,r]c) > 

-« E 

Oj ,0 2 CA;C 1 UC2=C 
O i naG 7 ;0;C i n(G»-) c #0 

x < (m dG - fhg G r ) ,d dGA \ dG K(- -> -;Ci\aG) d A \ aGtdG l Cl \dG conn .m(C 2 ; cr C2 , ??c 2 ) > ( 3 - 18 ) 
+ | ^ < rn(C 2 ;(Tc 2 ,ric 2 ), (naG^GrllaG) -1 rn(C 3 ;a C3 ,vc 3 ) > 



C 2 ,C 3 CA;C 2 UC 3 = C 
C i neG#0;O i n(G»') c #0 



9 
2 



E 



C 1 ,C 2 ,C 3 CA;C 1 UC 2 UC3 = C 
C i nSG 7 ;0;C i n(G'')<:^0 

x < rh(C2;o-c 2 ,Vc 2 ),d dGA \ aG TZ(- -► • ; Ci\dG) d A \ dGjdG l Cl \ dG conn .m(C 3 ; cr C3 , 77^) > 

with the short notation mgg = mg r 8 A,7 ' G "'' CTGr , The bounds are clear now from the bounds 
on the resolvent, the choice of r and the (trivial) control of the Cj-sums, i.e. provided by 

all subsets S 1 ,S 2 ,S 3 CC (3.19) 
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To proceed with the proof of Proposition 3.1 and high temperature-expand the .ff HT -terms 
we use the subtraction of bounds-trick to ensure the positivity of the resulting activities. We 
thus write for fixed G 

- ^ oca H Gr (mg a , <7 G r ,var;C,ac,vc) 

g Cn8G^0;Cn(G r ) c #« 

E„-a|C| 
CCAiCconn. to (CJj e 

II g CnSG 5 s»:Cn(G r )^« ^2 20) 

(G r )iconn. cp. of G r 

E cca (n(G r ,C)e- a W-H%g G r(.m aG ,a G r , Va r;C,a G , Vc )) 

X g CnaG#0;Cn(G r ) c 7i0 

where n(G r , C) is the number of connected components of G r that are connected to C (i.e. have 
(CT)i n C / 0). The exponential in the last line can then be cluster-expanded and gives 



Y\ cca (n(G r ,C)e-^ c ^-H^ ar (mg G ,a G r,T, G r;C,ac,Vc)) 

e cnaG^diCnfc)^* 

= X] PdG,GA m dG,^G-,VG-;K,a K ,VK) 

KCA;K = Or 

ifn8G#«,Kn(G r ) c #9 



(3.21) 



with < Pqq G r (mdG: cr G r : r lG r ')K,aK, , riK) < e "l-^l. Here we use the convention that 

PdG,G- {mdG,(^G-,VG-;K = Q,<t k ,Vk) = 1- 

Note that the resulting activities factorize over connected components of KUG r ; this is due 
to the (trivial) fact that the number n(G r , C) that enters the definition of the contour activities 
depends only on those components of G r that C is connected to. We put 

E„-°|C| 
CCA;Cconn. to (G r )i e 

p seo {dG,G') := [[ e cnea*« i cn(G'-) *« 

(G r )iconn. cp. of G r 

and note that 

l>p^(dG,G r )>e-\ Gr \ e ' c ° nst " (3.23) 

We can finally carry out the integral on dG to get the form as promised in the proposition. 
In doing so it is convenient to pull out a normalization constant and introduce the normalized 
Gaussian measures on dG corresponding to the Hamiltonian on the r.h.s. of (3.8), given by 



f (m 9gAG T,Og AG T),ri G r,a G r 

/ Hg,g- (dm dG )f(m dG ) 



_ Jdm dG e- AH ^ (maG) f(m aG ) 



.„(™ag A G'''»9 A G4 , 'C'»G'' / , n. 

f dm' dG e- AH 8oJr Kg) 
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So we can write 



dm dG e sg.a l m8G€ f/aeJ G 



= (27r) J ^ L Jdet(n8 G (o- g A G 0" 1 ne G ) P geo (9G,G 

" KCA;K = Or 

l 1 dG,G r {dm dG ) l mgG€U 9G p 9GjG r {m dG ,a G r,r] G r;K,a K ,r] K )I G 

(3.25) 

This has in fact the desired form (3.5) with the obvious definition 
/%a« (CG;^,^) :=p*°°(dG,G r ) 

x / 4S:^ Gr '° aAGr) ' , ' Gr ' CTGr (^G) l meGeC /ac pf G , G . (m^,^,^;^^,^)/™ 3 -™-*^ 

(3.26) 

with i"C = (5\G r on the r.h.s. Note that these activities factorize over connected components of 
G. 

In view of the trivial bound (3.23) on the geometric activity (3.22) and the normalization 
of the measure, the bounds follows from the HT-bounds and the bounds on the anharmonic 
activities I. The value of the 'Peierls constant' a is now clear from a = Const min{(2r + 
l)~ d log a}, assuming that both terms in the minimum are sufficiently large. 

To finish with the proof of Proposition 3.1 is now an easy matter. Using the formula for 
the determinant from Appendix (A. 3) we can write 

- ^ det (n aG (a - qA^y 1 n dG ) 

det(a-gA M eJ 

1 dPt (ll,, ( ;j„ <lA (; , ) ' II„ (l .) 



det ( 



o — ?Aa) det(n 9G (a- (? A A )- 1 n 9G ) 



x det (a - c/A G ) 



Remember that the correction given by the middle term on the r.h.s. stems from the lack of 
terms with range longer than r in the quadratic form of (3.24) that we had cut off. The random 
walk representation then gives the following expansion whose proof is given in the Appendix. 



Lemma 3.4: 



det (ll dG (a - qA Gr ) 1 Il aG ) _ 2 V CCA e det ( G ) 

/ _ g CnSG#»;Cn(G r ) c #9 (3.28) 



det (n 9G (a - qA A ) 1 n aG ) 



where < e det (C) < e a \ c \ with a ~ const log -. 
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Next we use subtraction of bounds as in (3.20) to write 

-E CCA £ det (C) 

g cnSG#e ; cn(G r ) c #9 

K<ZA;K = $ Or 

ffn8G#«,Kn(G r ) c #d 



xp^(dG,G r ) PaG, G r{K) 



KCA;K=I Or 
Kn8G#«,Jfn(G'')V« 



(3.29) 



where 1 > p geo - det > e -l Gr l e °°" s " and < /ogg iGr (A") < e -™" st °l^l. So we get 

li a i - |A| 1 _ ; n f , u^dA^A^A I i \ 

e b \ A \Z™ aA (a A ) = (2ir)*(det(a-qA A ))~e ™'a H * ™ 
x (2vr)"^ Vdet (a - qA G ) ^ p™^ (g, G; a G , V6 ) 

G:0CGCA G:Gca (3.30) 

G r CG 



This can be summed over G, G, K (collecting terms that give the same G U K) to yield the 
claims of Proposition 3.1. 



IV. The effective contour model: Gaussian case 

It is instructive to make explicit the result of our transformation to an effective Ising-contour 
model at first without the presence of anharmonic potentials where the proof is easy. In fact, as 
we will explain in Chapter V, the work done in Chapters II and III will then imply that a weak 
anharmonicity can be absorbed in essentially the same type of contour activities we encounter 
already in the purely Gaussian model. 

We remind the reader that in the purely Gaussian case the Ising-weights (t (^p A aA ' VA ^ (oa) 
are obtained by normalizing exp ^— inf mAeJR A h^a,va,ca ( mA )^ by its OA-sum. For simplicity 
we restrict now to the boundary condition fh x = m* for all x (that is everywhere in the minimum 
of the positive wells). 

We will now express the latter exponential as a sum over contour-weights. To do so we 
use the following (by now standard) definition of a signed contour model, including +-boundary 
conditions. 

Definition: A contour in A is a pair T = (T, o\) where TcA (the support ofT) and 
the spin- configuration a\ € { — 1,1} A are such that the extended configuration (o~a, +^z d \A) * s 
constant on connected components of Z d \T_. 

The connected components of a contour V are the contours I\ whose supports are the 
connected components T { of T and whose sign is determined by the requirement that it be the 
same as that ofT on F i . 
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A contour model representation for a probability measure v on the space {— 1, 1} A of Ising- 
spins in A is a probability measure N on the space of contours in A s.t. the marginal on the 
spin reproduces v, i.e. we have 



KK})= £ n({t}) 



r 



Recall that, in the simplest low-temperature contour model, arising from the standard 
nearest neighbor ferromagnetic Ising model, N ({r}) = Const x p(T) is proportional to a (non- 
negative) activity p(T) that factorizes over connected components of the contour and obeys a 
Peierls estimate of the form p(T) < e _T '-L There is a satisfying theory for the treatment of de- 
terministic models with additional volume terms for activities that are not necessarily symmetric 
under spin-flip, known as Pirogov-Sinai theory. For random models then, while the activities will 
be random, there have to be also additional random volume-contributions to N ({r}), even when 
the distribution of the disorder is symmetric, caused by local fluctuations in the free energies 
of the different states. The fluctuations of these volume terms are responsible for the fact that, 
even in situations where the disorder is 'irrelevant', not all contours carry exponentially small 
mass but the formation of some contours (depending on the specific realization) is favorable. It 
is the control of this phenomenon that poses the difficulties in the analysis of the stability of 
disordered contour models and necessitates RG (or possibly some related multiscale method). 

To write down the Peierls-type estimates to come for the present model we introduce the 
'naive contour-energy' (i.e. the d — 1-dimensional volume of the plaquettes separating plus- and 
minus-regions in Z d ) putting 

{x,y}cr,d(x,y)=l ^er.ygSA k • > 

d(x,y) — l 

again taking into the interaction with the positive boundary condition. 

Then the result of the transformation of the purely Gaussian continuous spin model to an 
effective Ising-contour model is given by the following 

Proposition 4.1: Suppose that q is sufficiently small, q{m*) 2 sufficiently large, a is of the 
order 1 and 5 < Const m* with a constant of the order 1. 

Then there is a a ^-independent constant K\ (tja) s.t. we have the representation 

xe Eccv +( . A) ^ Ufi (-)-E ccv - ( „ A) ^(,c) £ (4.3) 

r 

°A( r )=< 7 A 
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for any a a, with V ± (o-^) = {x G A; o~ x = ±1}. Here 

(i) rjc ► S^ uB (rjc) are functions of the random fields indexed by the connected sets CcA. 
T/iey are symmetric, i.e. S^, auB (—r]c) = — 5g auB (rjc) and invariant under lattice- shifts. For 
C = {x} we have in particular S^ uB (rj x ) = ^"yidq rjx ■ 

(ii) The activity po(T;nr) is non-negative. It factorizes over the connected components ofT, 
i.e. 

Po(T;7]r)= ]J PoFiiVTi) ^ 

Fi conn cp. of r 

For T not touching the boundary (i.e. c^aE = %) the value of po(F;rjr) is independent of A. We 
then have the 'infinite volume properties' of 

(a) Spin- flip symmetry, i.e. p ((F, cr A ); rjr) = Po((L, -o"a); -Vt) 

(b ) Invariance under joint lattice shifts of spins and random fields 
Peierls-type bounds: There exist positive constants /3 G auB>/3 s.t. we have the bounds 

< p (r;r?r) < e -^( r )-^a uE |r| (4 . 5) 
uniformly in r]r_ S [—5,(5]— where the 'Peierls- constants' can be chosen like 

^ ^(^W ^^C^xminjlogi,,™-^) J - (4.6) 

The non-local random fields obey the estimate 

\S2T u& (vc)\ < 5m*e-^ c \ (4.7) 
for all \C\ > 1 with (3 = Const x log -. 

Remark 1: This structure will be familiar to the reader familiar with [BK1] or [BoK] (see 
page 457). Indeed, the above model falls in the class of contour models given in (5.1) of [BK1] 
(as written therein for the partition function). This form was then shown to be of sufficient 
generality to describe the contour models arising from the random field Ising model under any 
iteration of the contour-RG that was constructed in [BK1]. (The additional non-local interaction 
VF(r) encountered in [BK1] is not necessary and could be expanded by subtraction-of-bounds 
as in (3.20), giving rise to enlarged supports T, as it was done in [BoK]). 

Remark 2: There is some freedom in the precise formulation of contours and contour activities, 
resp. the question of keeping information additional to the support and the spins on the contours. 
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[BK1] speak of inner and outer supports, while in [BoK] it was preferred to define contours with 
activities containing interactions. The latter is motivated by the limit of the temperature going to 
zero (making the interactions vanish). Since we do not perform such a limit here, we present the 
simplest possible choice and do not make such distinctions here, simply collecting all interactions 
from different sources into 'the support'. 

Remark 3: The magnitude of [5 ~ Const qm* 2 is easily understood since it gives the true 
order of magnitude of the minimal energetic contribution to the original Hamiltonian of a nearest 
neighbor pair of continuous spins sitting in potential wells with opposite signs. This term appears 
again in the estimate on /3 GauB (up to logarithmic corrections) together with a contribution of the 
same form as (3$. The latter comes from a straight-forward expansion of long-range contributions. 
The last term in (4.6), m*S, is a trivial control on the worst realization of the random fields; it 
could easily be avoided by the introduction of so-called 'bad regions'. These are regions of space 
where the realizations of the random fields are exceptionally (and dangerously) large in some 
sense and, while comparing with [BK1] or [BoK], the reader might have already missed them. 
Indeed, a renormalization of the present model will immediately produce such bad regions in 
the next steps. Of course, we could have started, here and also in the presence of anharmonicity, 
with an unbounded distribution of the r/ x . In the latter case we would have to single out 
regions of space where the behavior of our transformation to the Ising-model gets exceptional 
(i.e. because we lose Lemma 2.5.) We chose however not to treat this case here in order to keep 
the technicalities down. 

Proof: An elementary computation yields the important fact that the minimum of the quadratic 
Hamiltonian (2.6) with any boundary condition m is given by 

- inf ^»a,^ A( ) = _^2 < fl + £(^£| A | 
m A eiR A 2q 2 



am 



<r] + rj d(A c ) (qm),R A a A > A 

1 q 

- 2~ < V + Vd(A°)(qrh),R A (rj + fjd(A*)(qm)) > A +- ^ 

" x£A;y£dA 

d(x,y) = l 

with f}Q( A c^(m) := d At dA?ndA denoting the field created by the boundary condition. We subtract 
a term that is constant for a A (and thus of no interest) and write 

inf tf™eA,*7A,a A (mA) _ inf H m 9A , VA ,l A (mA) _ « < VAiRa i a >a 

m A m A q 

a 2 (m*) 2 . , am* n (4.9) 

= ~ (< cta, Rao~ a >a - < 1a, -KaIa >a) < Va,Rao-a >a 

2q q 

- am* < fj 9 ( A c)(m),RA (o~ A - l A ) > A 

The first term on the r.h.s. gives rise to the low-temp. Peierls constant; the next term is a 

weakly nonlocal random field term (suppressed by the decay of the resolvent) and the last term 
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the symmetry-breaking coupling to the boundary. As in Chapter III we use the random walk 
representation R A = ^ CcA Ti (• — ► ■ ; C) (see Appendix (A.ll))and decompose according to the 
size of C's. 

As the first step for the contour representation we associate to any spin-configuration a\ G 
{ — 1,1} A a preliminary (or 'inner') support in the following way. Choose some finite integer 
r > 1, to be determined below, and put 

£a(°a) := {x € A; 3y G A s.t. d(x,y) < r where a x / <7 y } 

(4.10) 

U {x G A; <i(x, <9A) < r + 1 where = — 1} 

The second term makes this definition A-dependent by taking into account the interaction with 
the boundary leading to the (desired) symmetry breaking for contours touching the boundary. 
For given a a the activities po(T; r]r) to be defined will be non-zero only for supports rDr + (<7A). 
The range r will be chosen below in such a way that the terms corresponding to interactions with 
range larger than r have decayed sufficiently so that they can be high-temperature expanded in 
a straightforward way. This choice then also determines the value of the Peierls-constant for the 
low-temperature contributions. 

Keeping the small C's of diameter up to r define the (preliminary) 'low-temperature activ- 
ities' 

p LT '™ 8A (<7 A ) 



Zc 



q- i (m*)^ 



{<crc,K{- -►■ ;C)a c >-<lc,n{- -►■ ;C)lo>)+am*<tja(A=)(™),R(' -►■ ;C)(a c -l C )> 

(4.11) 

Note that the 'inner support' (4.10) can be trivially rewritten as 

E + K) = U cu U o 

CCA;diam(C)<r C conn, to SA v ' 

< T 0^ 1 and a C^- x C diam(C)<r;CT c #l c 

which shows that it is just the union of all connected C's with diameter less or equal r that give 
any contribution to the sum occuring in the exponent of (4.11). So we can rewrite 

e -inf mA H™ 9A ^ A '" A (m A )+inf mA H-8A^A,iA( mA ) + £sl <r;AijRA i A>A 

^TT" CCA;diam(C)<r <VC — ■ ;C*)<T C > 

_ p LT,msA (cjA)e »c^ c «" st 

e ^f Eccv +( r)<^' K (- - ;C)lc>-^ Eccv-(r)<^' K (- ^ ^)^> (4 



° 2( r g '' )2 (< CT c,R(- ;C)a c >-<lc,n(- ;C)l c >) + ii f :1 <r;c,K(- ;C>c> 



^ CCA;diam(C)>r 
g (Tfj^const 

CCA;diam(C)>r «™* < »?S( A<= ) (™) ' 5 C) (<T C - 1 C ) > 

e cna(A) c #o 
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The terms in the first line depend only on quantities on T + (crA) and factorize over its connected 
components. They will give contributions to the activities p$. The terms in the second line are 
the small-field contributions to the vacua given by 

S^ B (vc):=^-<Vc,n(- -■;C)lc> (4.14) 

The terms in the last two lines are small (since only C's with sufficiently large diameter con- 
tribute) and only non-zero for C's intersecting with r + (cTA) or touching the boundary. They 
can be expanded. 

Let us see now what explicit bounds we get on the low-temperature activity (4.11). Keeping 
only C's made of two nearest neighbors x, y = x + e we have the upper bound 

£ ^j-^(<crc,n(- ^.;C)a c >-<lc,K(- --;C)ic>) 

CcA;dia m (C)<r Q 

a 2 (m*) 2 ( ^ 

<- ° [m ' £ U{x ^ y-C = {x,y})l ax ^ y 

q {x,y}cr+(cr A ),d(x,y) = l 

Computing 

1 °° / 1 \ k 1 

U (x -+ x + e; C = {x, x + e}) = _ ^£ j = (c + 2d)2 _ 1 (4.16) 

with c = a/g we get an upper bound on the l.h.s. of (4.15) of -2(3 Y,{x, y }cr+(cr A ),d(x,y)=i l <y^ y 
where (3 is given by (4.6). Applying a similar reasoning on the boundary term, thereby using 
that 1Z (x — ► ; C = {x}) = c _^ 2d , gives the bound 

£ am* < Tj a{ Ac)(m),K(- -> • ; C) (a c - lc) > 

CcA;diam(C)<r 

„ 2 2a ^ (4.17) 

< — aim ) > lo- __1 

~ yv a + 2dq ^ 

d(a;,y)=i 

Since the modulus of the prefactor in the last line is larger than 2/3 we get an energetic suppression 
of 

^T.maA^^j < e -2f3E s {r+(a A ),a A ) < e -/3£ s (r+ ( ( 7 A ), CTA )-/3(2r+l)- d |r+ (<x A ) | ( 418 ) 

Using Yl yR\;x, y < 1/c for the next term in (4.13) we have immediately 



am 



£ <r/c,^(- ^-;C>c ><m*5|r+(a A ) 



9 — (4.19) 

^ CCA;diam(C)<r 

This finishes the Peierls estimate for the low-temperature contributions. 
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Let us come to the treatment of the 'high-temperature parts' in (4.13) now, proceeding 
algebraically at first. Using subtraction-of-bounds as in Chapter III (3.20) we get the high- 
temperature expansion 



CCA;diam(C)>r 
g (TQ^const 

= /5 ge °(r + (a A )) 



- TTp («rc,K{- ;C)a c >-<lc,'R.(- ;C)l c >) + ^<r; c ,K(- ;C)ct c > 



E 



p HT1 (K,a K ,r, K ) 



(4.20) 



K"CA;diam(K)>r OI ft" - 



if the terms in the exponential on the l.h.s. are sufficiently small. To control them we just use 
the bound (A. 12) 

1 / 2d \ |c|_1 



n{x^ y -c)< 

yez d 



This gives the deterministic bound upper bound on the first two terms in (4.13) of 

a 2 (m*) 2 



2q 



if we have 



<a c ,n(- -^■;C)a c >-<lc,K(- • ; C) l c >| < e" Q|c| 

log(am* 2 ) 



(4.22) 



1 - 



(|C|-l)log(l + 4 



(4.23) 



which is in turn bounded by a := | log ^1 + — \ for the C's in the above sum if we put 



; log (am* ) 



+ 1 ~ 4 



logm* 



(4.24) 



Remember here that we are interested in the regime of - small and m* 2 even larger. 



Assuming (4.24), with < S the random field contribution is estimated by 

-e 



^\< VC ,K(. ^.;C)a c >\<^-^ c \ 



(4.25) 



q am* 
where can use that < Const. The estimates on 5 , g aufi (r/c) are obtained in the very same 
way. 



In passing we verify that all activities constructed so far are invariant under joint flips of 
spins and random fields (inside A). The boundary terms can be expanded similarly giving 

ccA;diam(c)>r am" <f] d{ A c ) (m) ,K( ■ ->• ;C) {a c - 1 c)> 

g Cna(A) c ^« 

-2 X] CcA;diam(C)>r am* <77 0( A c ) (m) ,1Z( ■ — >■ ;C)lc> ^ ccA;diam(C)>r am* <?7 S ( A c ) (m) • — >■ ;C)(ct c + 1 c )> 
= g Cna(A) c 7i0 g cnO(A) c 7 i0 

-2 X] ccA;diam(c)>r am* <ij a{ A c ) (m) ,K{ ■ ;C)1 C > 
= g Cn8(A) = #« 



E 



KCA;diam(K")>j- Or K = 

/fn9(A)«^« 



(4.26) 
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This gives 

-inf mA ^ aA, ' 7A '' JA (m A )+inf„ lA H™ sa ' v *' 1a (m A )+ ^<»? A ,.R A l A > A + 2 £ cc A ; diam(C) «m*<»j 8 (Ac)(m),R(^>;C)lc> 
e cna(A c )^0 

E CCA;diam(C)<r <VC W ;C)<T C > 

— p LT ' msA (cr A )e »c^ co " sl 

xf°(r> A )) ^ p HT1 (^,a x ,^) 

XCA;diam(K)>r OI K = 

K 1 CA;diam(K 1 )>r Or K 1 = 
K 1 n3(A)<: 7 £0 

(4.27) 

which proves the desired representation (4.3) with the obvious definition 
Po(r;f?r) 

^77~ ^ CCA;diam(C)<r <1C.R( I ~ > ' ;C)<7C> 

x r ° (r+ (a A )) 2 -° HT1 ( K ° ' a ^ ' ^° ) p HT2 ^ ' ' ^ ) 

K ,K 1 CA;K UK 1 ur+(u A )=r;diam(K i )>r Or K i= 
o K() ^consi ,K 1 nS(A) c ^0 

(4.28) 

The form (4.6) of the Peierls constant (3 GauB is now clear from j3 GauB = Const min{/?(2r + 
l)~ d , a } — m*5, assuming that both terms in the minimum are sufficiently large to control the 
entropy in (4.28) and the slight modification in the exponential bounds on p HT1 arising from the 
subtraction of bounds. 

V. The final contour model - Proof of phase transition 

We put together the results of Chapter III and IV to obtain the contour representation of 
the full model. It is of the same form as the Gaussian model of Chapter IV, while a modifaction 
of the Peierls constant (3 accounts for the anharmonic contributions. More precisely we have 

Proposition 5.1: Assume that the anharmonic I -weights (2.17) satisfy the Positivity (2.19) 
and the uniform Peierls Condition (2.20) with a constant e. Suppose that e is sufficiently small, 
q is sufficiently small, a is of the order one, q(m*) 2 sufficiently large. Suppose that 5 < Const m* 
and \U\ < Const m* with constants that are sufficiently small. 

Then the measures T ^/x^ m laA ' VA ^j on { — 1, 1} A have the contour representation 
T (MA m * l8A '" A ) K) 

contour , A r 

»A( r )=°"A 
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with the contour-model partition function 

*&ZL,A = E eEccv+(rA)5c(r,c) " Eocv - (rA,Sc( ' ?c V(r; ?? r) (5.2) 

r 

For the partition function (1-4) we have Z^ m laA ' VA = C A ~' VA Z^~ o '^ ui A with a trivial constant 
containing the contributions of Gaussian fluctuations that satisfies, a.s. 



lim i-logC+' ??A 
A~[Z d |A| A 



1 1 



(5.3) 



2 

The quantities appearing in (5.1) are as follows. 



(a-qA zd ) 1 --[log(a-gA zt i)] 0j0 -6 + -log(27r) 

- 0,0 — ■ 



(i) rjc i— ► Sc(r)c) are functions of the random fields indexed by the connected sets CcA that 
are symmetric, i.e. Sc{—f]c) = —Sc(rjc)- In particular we have S x (rj x ) = ^J^ d( J ]x- They 
obey the uniform bound 

\Sc(vc)\<m*5e- a f'^ c] (5.4) 



for all C with a fina i = const min jlog ±, log \ (J^*) } ' 



(ii) The activity p lBins (T;r]r) is non-negative and depends only on the indicated arguments. It 
factorizes over the connected components (as in (4-4))- ForT_ not touching the boundary it 
does not depend on A and has the infinite volume symmetries of a) invariance under joint 
flips of spins and random fields and b ) invariance under lattice shifts. 

There exist (large) positive constants (3,f3 s.t. we have the Peierls-type bounds: 

p Isi „ g (r ;7 7r)<e-^( r )-fel (5.5) 

uniformly in rjQ. Here [3 = q(m 2 ( a +2d q ) 2 - q 2 * s ^ e same as ^ n (4-8) and 



(3 = Const x min < 



1 2 / log - \ d I I log - 

log -,qm* — , log- — I V — m*5 (5.6) 

q \ log m* I e \ log m* ' ' 



Proof: Assuming the control of the anharmonicity, summarized in Proposition 3.1, the proof 
is easy. For any fixed a a we can cluster-expand the last sum in (3.1). Dropping now the 
dependence on the boundary condition rhdA = +m*lgA hi the notation we have 

log Yl P(Q o-g,Vg)= E e(C;a c ,Vc) 

G:0CGCA C:0CCcA 

= HC;ic,vc)+ £ HC;-ic,vc)+ £ HC;<jcvc) (5-7) 

CCV+((J A ) CcV-(cr A ) C(ZA;a c ^const 
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where the sum is over connected sets C and we have the bounds |e (C; ac, Vc)\ < e -consta\c\ 
with a given in Proposition 3.1. Together with the representation (4.3) for the purely Gaussian 
model this gives 

-inf c1r a iC 0A ' T,A ' <TA ( m A ) -fsi \ t \ 

e -A eHA 2^ p(G;a G ,VG) = K A {r] A ) 

G:0CGCA 

xe Eccv +( „ A )( S § aU ^)+^^^ (5 . 8) 

r 

cr A( r ) = ' T A 

Note that the C's in the exponential in the last line are in particular connected to T. Using 
subtraction-of-bounds as before we can expand those terms and, as we did before in Chapter III 
and IV, rewrite the last line in terms of a new (and final) contour summation as 

e Ec^». t ^) £ po^r) = £ p(T ;7 r) (5 . 9) 

r r 

o-A< r ) = < : 'A °"A( r ) = IJ A 

The values of the Peierls constants for the final activities on the r.h.s. follow from the statements 
of the Propositions 3.1 and 4.1 with a slight loss due to the control of entropy. 

Finally, to see the statement for the free energy, we start from (3.1) and recall the con- 
struction of the activities in the purely Gaussian case, starting from (4.16). Using the explicit 
expression (4.11) for the energy minimum in the Gaussian model in terms of the resolvent we 
obtain, with some trivial control on boundary terms, using the SLLN applied on the random 
fields the desired formula 

- ~ Sfc ff *"'"" ,,A <"»»> " i{S. m logdet(o ~ ,Aa> ]og{2 ' ] (5 - 10> 

IEti^ r i 1 1 

= if- [(a ~ q^z-y'l Q Q " g [ lo S ( a " 9 A *««)]o,o ~ b +2 log(27r) 

<> 

The following result provides control of the original measure in terms of the coarse-grained 
one up to two corrections: 

Proposition 5.2: Assume the conditions of Proposition 5.1 and suppose that ttiqa G (U + ) dA . 
Then we have 



,.m a A,VA 



m 

~ 2 



< (t (a*a oa '" a )) k«o = -i] + e - constQ + e - cons * (m * )2 



(5.11) 
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where a = const x min |log |, log \ ( log^* ) | * s 9^ ven ^ n Proposition 3.1. 

Remark: The first term on the r.h.s. accounts for the anharmonicity, the next one for the 
Gaussian fluctuations. 

Proof: We carry out the transformation that led to Lemma 2.3 while carrying through the 
indicator function 1 ^ m * to get 

I dm A l m <m , e -<--^-) = e -|AIV e - inf ^^ A ^ A, " A ^ 

JffiA m -0<— ^ 

o A 

£ (2^ (d«(.-,A AXC ))-* / * B , e-^— 

!:GCA •* 

+ ^ (2vr)^ (det(a-,A A ^)) _i ! dm dG e- AH ^X "^1^^^/, 



T m gA ,m gG ,r]G,o-G 



G 



G:GCA 
8G3i 



G_GCA 



x ji~n gA ,m0G,r]G,VG 
G 



(5.12) 



with /g;x = ^g-L ~~ ^G-L ( su P erscr ipts are dropped now) where we have defined 



(5.13) 



2 



We use the same notations without the subscript xo on the l.h.s. to denote the integrals without 
the l m <2 : on the r.h.s. so that we have Iq = Iq — Iq- Note that it is not clear anymore 
that Ig;x is positive for any sign a Xo and dominated by Iq- To bypass this little inconvenience 
we argue as follows. Let us slightly enlarge b in Chapter II by putting a factor 2 in front of 
the fraction of integrals in the definition (2.42). This leaves b very small and all subsequent 
arguments based on a fixed choice of b remain valid. Going back through Lemma 2.4. we see 
that this definition implies that even 

7 (2) < 2 -l G l/ G 1) (which can be seen as a strengthening of 
the positivity of Iq). But from this we have in particular that 

lG; X0 = /gLo " ^gL * J gL <^<Vg (5-14) 

We use this estimate on the last G-sum in (5.12) and bound the second G-sum in (5.12) by the 
corresponding expression without the indicator. Carrying out the mac-integral as described in 
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Chapter III we get from this the bound 



G:GCA •> 
G3x 

+ £ (2.)^(det(a- g A A ^)p/^ 



dA,m aG ,ri G ,a G 
x 



GGCA 
8G3i 



i Tm aA ,m gG ,ri a ,o- G 
m eG eU 9G 1 m :cn <^ 1 G 



< 2- (2vr)- (det(a-gA A )) 1 ^% G (G;ff G ,r, G ; 

G:x eGcA 



(5.15) 



Using the positivity of the activities in the last line we can use the usual Peierls argument on 
the fixed-cr contour model appearing in (3.1) that controls the anharmonicity. So we estimate 



]T ffi l °e\ G (G\aG,f]G) < £ P* 88A °(Go;a Go ,f^ ) £ ir B °* a (G;* G ,ri G ) 

G:io£GcA G„:i„eG7cA G:GcA 

G:GcA 

(5.16) 

where the first sum is over connected sets Go and we have used Proposition (3.1) for its estima- 
tion. 

To treat the first G-sum in (5.12) we note that the expectation outside the anharmonic 
contours is given by the one-dimensional Gaussian probability, 



dm dG e- AH oSr A ^o h < ^ = (27r) ^(detfa-,A / ^ 7 



*A\G 



X M 



rnoA,m aG , VA .-,<T A .- 



m 



A\dG 



-1 

x ,x 



(5.17) 



\ m < m * ) 
v I'I-xq i 2 



with the notation J\f[a; <J 2 ]{4>) = f_ 



CO e la' 2 - <j>(x) 



. We use the uniform control on the expecta- 



-OO V27TO- 2 

tion value given by Lemma 2.5 and the fact that the variance occuring in (5.17) is of the order 
one, in any volume. If a xo = +1 we have from this, uniformly in all involved quantities that 



■m gA ,m aG ,r] . -^er — 



m 



A\dG 



-1 

Xq,Xq 



(U<i^) < e 



-const (m* ) 2 



(5.18) 
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so it can be pulled out of the mac-integral. For a Xo = — 1 we use the trivial bound 1 to write 



dm dG e G,A {maG) l maG€U ec / dm A ^e a\o I 1 7 



&GCA 



T-m gA ,m0G,riG,<7G 



< (e- co " st( ^ )2 l^=i + l CTiC0= -i) 

GiGCA ^ 

< (e- C ° nSt(m * )2 +l CT;co= -l) 

G:GcA ^ 

(5.19) 

Now it is simple to put together (5.12), (5.16)-(5.19) and rerunning the next steps of the trans- 
formation yields the claim. 

Applying the information of [BK1] we obtain the main result of the paper. 

Proof of Theorem 1: We apply statement Theorem 2.1 [BK1] on the measure T (jj% m * laA ' r ' A ^ . 
Indeed, this is justified from Proposition 5.1 which implies that this measure is contained in the 
class of contour measures described in [BK1] Chapter 5 'Flow of the RGT', Paragraph 5.1. We 
note that of the three constants (3, /3, a fi na i (controlling the exponential decay of the activities 
in terms of the volume resp. in terms of the naive contour energy, and the decay of the non-local 
fields) the constant (3 is the smallest. 

So statement (2.3) from [BK1] gives in our case that for d > 3, (3 large enough and a 2 small 
enough we have that 



IP 



T ^X m ' ldA ' nA ) Wx = -1] > e- const ^ < (5.20) 
We apply our Proposition (5.2) and note that the two correction terms given therein are also con- 
trolled by e ~ const 13 (with possible modification of const .) From this in particular the estimates 
of Theorem 1 follow. 



Remark: We have not given an estimate on the value of 7 as a function of q and m* . This 
would of course follow from a more careful estimate of the best value of the 'anharmonicity- 
constant' e (which is entering (3) as a function of q and m* (see Chapter II) and is left to the 
reader. 



Finally, Theorem 2 for the -theory follows immediately from 
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Proposition 5.3: Assume that the anharmonic I -weights (2.17) satisfy the Positivity (2.19) 
and the uniform Peierls Condition (2.20) with a constant e (that is sufficiently small). Let 
be any continuous spin Gibbs-measure obtained as a weak limit of [j,™ eA ' VA along a sequence of 
cubes A for some (not necessarily positive) continuous- spin boundary condition rh G U z . 

Then the measure T(/i^ c ) on {— 1, l} zd is a Gibbs measure for the absolutely summable 
Ising-Hamiltonian 

m sm9 M 

a 2 (m*) 2 v^. » — i » \ — l ^ , \ (5-21) 

= 2 Z^^ a ~ qA ^x, y (J xO-y-am 2^(a-qk zd ) x y n x (7 x + ^ $c (pcWc) 

x,y x C:\C\>2 

where the interaction potentials <&c(o~c,Vc) = ^c{—o~c, —"He) obey the uniform bound \<&c{o~c-> Vc)\ < 
e -const a \c\ j or dl c with a = CQnst x min | log i >log I (Jg-L-) | as in (3.3). 

Remark: Note that it follows in particular that the interaction will be the same e.g. also in 
continuous spin Dobrushin-states [Do2] (that are believed to exist) one could construct using 
the boundary condition +m* in the upper half-space and — m* in the lower half-space. 

Proof: Denote by v (ay) the usual restriction of (5.21) to the finite volume V, obtained 

by keeping the sums over sets {x,y} and C that intersect V and putting the spin equal to a Z d 
for x $ V. Following [BKL] it suffices to show that, for each <f z ,i we have that 

Km lim '(^W) = e ^ 

along (say) sequences of cubes where 



JlRAl ,6A 2 



(5.23) 



This is clear, since (according to our assumption of weak convergence) there is a subsequence of 
cubes Ai s.t. the inner limit exists and equals (T(/i^ c )) (cry |<7A 2 \y) . Summing Proposition 3.1 
over the spins in Ai\A 2 we have then 

zZ^(a A2 )=e- b ^(2n) L ^(det(a- q A Al ))- 1 * £ e"^^ 

<5- Al \A 2 

x Z^j p m9aAlG {G;a GnA2 ,d- GnAl \ A2 ,ri G ) 

G:0CGcAi 

(5.24) 

From here the proof is easy, given the explicit formula (4.11) for the minimum and the absolute 
summability of the polymer weights, uniformly in the spins and random fields. 
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For the convenience of the reader we give a complete proof for the simplest case of vanishing 
anharmonicity w(m x ) = 0, and vanishing magnetic fields rj x = 0; it illustrates the way boundary 
terms are entering. Using (4.11) we have indeed 

!= Const X ^ e ^-^^<(^A 2 ,cr AlNA2 ),itA 1 (CTA 2 ,cr AlXA J> Al +am*< } 7g (Al c ) (m),i? Al ( CTA2 , C r AlXA2 )>A 1 

(5.25) 

Now, using the exponential decay of the resolvent, 

„m eA 1 ,r]A 1 f - \ 
Z A 1 {<7V,0-A 2 \V) 

— Const X e ° ^ 1 <{ a v ,°a 2 \v),Ra 1 (crv,^A 2 \v)>A 1 +am* <fjg (Al c ) (m),R Al (w ,<?a 2 \v)> A± 

x g ° ( ™ ) <(<?V ,0"A 2 \v)i-Ra 1 0-A 1 \A 2 >A! + ° ( ^ ) <g"A 1 \A 2 ,-RAj^ <T Al \A 2 >A X +am* < J7g ( Al c ) (ttOjHAj ^A^ \ A 2 >A 1 

0-A!\A 2 

= Const X e ^ I i ni <K^A 2 \v),«A 1 (av,SA 2 \v)>A 1 ±Con S t|A 2 |e-"' d ' st ( A ^ A ^ y e ±Con S t | V | e — 'd»t(V,A|) 

x 'y * g ° ( ™ - <P'A 2 \V;-Ra 1 'TA 1 \A 2 >Aj g+ ° ^ - < cr A 1 \A 2 i^A^A^Aj > A! +am* < rjg( Al c ) (m) ,Ra 1 f Al \ A 2 >Aj 
C Al \A 2 

(5.26) 

The terms in the last sum do not depend on ay so that we get 

Z A 1 {<7V,VA 2 \v) 



EryKldA, ,VA 1 (~ \ 
~a v Z A, {VV,0-A 2 \V) 

, n ... , -«'dist(A 2 ,A9), „ -«'dist(V,A§) g" ^ 1 <(av,VA 2 \v),RA 1 {<?V,VA 2 \v)>A 1 

_ g±Const \A 2 \e ^ i'±Const \ V \e v 2' 



V p " ^ ) < (^v-,o-A 2 \v),flA 1 (ctv,cta 2 \v)>a 1 

(5.27) 

with uniform constants. Taking first Ai \ /Z d (using that R\ 1 \ A ^^ Rz d \^ 2 ) an d then A2 f 
we get in fact the desired result in our special case. 



The (random) non-Gaussian case follows easily from the cluster expansion of the G-sum 
in (3.1). Indeed, since we have uniform exponential decay of the activities ^" 8 sa g , the cluster 
expansion gives us quantities <&™° oaC (ac'iVc) that obey a uniform bound of the form as desired 
s.t. we have 

p mo aAO(G;a G ,^ G ) = e^conn,0cccA*c ^o;vc) 

G:0CGcA 

The Gibbs potential in (5.21) is then given by the value of $c f° r polymers C that are not 
touching the boundary. With estimates on boundary terms as in (5.26) the claim (5.22) follows. 
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Appendix 

For easy reference we collect some formulae about quadratic forms and the random walk 
expansion of determinants and correlation functions we use. We start with 

Lemma A.J: Let Q A be symmetric and positive definite. Let VcA and write, with obvious 
notations, 

Qv Qv,A\V 



Qa\v,v Qa\v,a\v 



Then we have the following formulae, 
(i) 

Qt = 

[Qv — Qv,a\vQ a\v,a\v^ *a\v,v*J —Qv 1 (Qa\v,a\v — Qa\v,vQv 1 Qv,a\v) 

~Qa\va\v (Qv ~ Qv,a\vQ A \v a\vQ a \v,v) (Qa\v,a\v ~ Qa\v,vQv 1 Qv,a\v) 

(A.2) 

(ii) 

det Q A = det (UvQ^Uy) _1 x det Q A \ V (A.3) 

(Hi) For any z A we can write 

i < m A , Qaw-a >a - < m A , z A > A 

= X - < (m v - m z A A | J , (UvQ^Uy) _1 (m v - m A A | J > A (A.4) 

2 \ mA - V ~ m A\V ) ' <5a\V - m A\V J >A~^< Z A , Q A Z A > A 

where the 'global minimizer' m A A = Q A 1 z A is the minimizer of the total energy, i.e. 

m A !-»■ ^ < m A , Qa^a >a - ^ < m A, z a >a ( A - 5 ) 
We write Qv = ^vQa^v , Qa\v,v = ^a\vQa^v ■ The 'conditional minimizer' 

m A A \v' mV = Qa\v ( Za \ v + QA\v,vm v ) (A.6) 
is the minimizer of the function 

m A \ v ^ ^ < ( m A\v,™v) ,<5a (m A \ v ,m y ) > A -- < (m A \y,my) , z A > A (A.7) 
for fixed mv ■ 
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Remark: The quadratic forms on the diagonal of the r.h.s. of (A. 2) are automatically positive 
definite. 



The proofs are easy and well known computations and will not be given here. Next we 
collect some formulae and introduce notation concerning the random walk representation. 

Lemma A. 2: Denote by 1Z the (non-normalized) measure on the set of all finite paths on Z d 
(with all possible lengths), defined by 

i \ ItI+i 

J_) (A.8) 
c + 2d) 

for a nearest neighbor path 7 of finite length \^\. Then we have 

R\-,x,y = Tt{l from x to y; Range{^)<Zh) (A. 9) 

where Ranged) = {7^ t = 0, . . . , k} is set of sites visited by a path 7 = ("ft)t=o,...,k of length 

h\ = k. 

Proof: Write A® = 2d — T A where T A . X:V = 1 iff x, y € A are nearest neighbors and T A;Xty = 
otherwise. Then 

00 / \ \ t+1 

R A , x , y = (c + 2d-T A )^ y = £ [ 7T ^) (tf) Xiy (A.10) 

which proves (A. 10). 

We will also use the obvious matrix notation 

(K (■ -> • ; C)) x y = n (7 from x to y; Range( 7 ) = C) (A.ll) 

so that one has the matrix equality Ry = J2c<zv ^ (' — * ' > C) for any volume V. We need to 
use a bound on its matrix elements at several places. Let us note the simple estimate 

1 / 2d \ I ' -1 

K (x -> V ; C) < U (7 starting at x, length( 7 ) > |C| - 1) = - I — ^ J (A. 12) 

We will use these notations at many different places. As an example, let us prove formula (3.11). 
Indeed, we have 

(n 9G (c - A A ) 1 n 9G ^ = c - (A dG + d dGiA \ dG R A \ dG d A \ dG:dG ) 

= C — (Ag G + dd G ^ A \9 G R G r\ dG d A \d Gt 9 G ) — dg Gt A\dG {Ra\8G ~ RG r \do) 9 A \dG,dG 

= (j^dG (c - A G r) 1 U dG ^j - d dGiA \ dG (R A \ dG - R G r\ dG ) d A \ dGidG (A. 13) 

= (jldG (c - A G r) 1 Uqg^ - d dGjA \ dG TZ(- — > • ; C) d A \ dGjdG 

OCA\8G 
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Here we have used Lemma A.l(i) in the first and third equality and Lemma A. 2 in the last one. 
Finally we give the 



Proof of Lemma 3.4: The random walk representation of the determinant is obtained 
writing 

log det (c — Ay) = log det 



= I V\ log (c + 2d) + Tr log (l + —^-jT v 

\ c + 2d 

and expanding the logarithm. Using (A. 3) we can then write 



(A.14) 



det (U dG (a - qA G r) 1 n aG ) det G ._ aG (c - A G ._ aG ) det A (c - A A ) 

log t f = log 



det(n 9G (a-gA A )- 1 n 9G ) " det G -(c-A G r) det A _ aG (c - A A _ aG ) 
°° 1 1 

= S J ( c + 2d)* ( TrG " (Tcr) * " TrGr " 9G C^-ac)* - Tr A (T A )* + Tr A _ aG (T A _ aG )') 



(A.15) 



t 

It is not difficult to convince oneself that we have that 

\t m Irn \t m. /m \t , m /m \i 



Tr G r(T G r) — Tr G r_ aG (T G »-_ aG ) -Tr A (T A ) + Tr A _ aG (T A _ aG ) 
= - #{7 : x ^ x; Range(7)cA; Range(7) n dG + 0; Range(7) n A\G r + 0; | 7 | = t} ( A ' 16 ) 



x€A 

So we get the form (3.28) putting 



1 1 

2e d -(C) := J] --— — J2 #{7 = x ~ x; Range( 7 ) = C; | 7 | = *} (A.17) 

t=2 ^ ' xec 

From this the bounds of the form e det (C) < e ~ const ( lo s ^ICI are c ] earj assuming that c is large. 
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